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ABSTRACT 

Using the Einstein-Hilbert approximation of asymptotically safe quantum gravity we present 
a consistent renormalization group based framework for the inclusion of quantum gravi- 
tational effects into the cosmological field equations. Relating the renormalization group 
scale to cosmological time via a dynamical cutoff identification this framework applies to 
all stages of the cosmological evolution. The very early universe is found to contain a pe- 
riod of "oscillatory inflation" with an infinite sequence of time intervals during which the 
expansion alternates between acceleration and deceleration. For asymptotically late times 
we identify a mechanism which prevents the universe from leaving the domain of validity 
of the Einstein-Hilbert approximation and obtain a classical de Sitter era. 



1 Introduction 



It is generally believed that in the very early universe, at times smaller than the Planck 
time, quantum gravity effects play a crucial role. Since the physics in the Planck era 
prepares the initial conditions for the subsequent classical evolution of the universe and, 
as a result, is likely to be responsible for many of its features we observe today it is clearly 
very desirable to gain some understanding of the quantum gravitational processes which 
took place immediately after the Big Bang. 

The construction of a consistent quantum theory of gravity which is indispensable for 
investigations of this kind is a major challenge, of course. In trying to set up such a 
theory the first question to ask is what are the degrees of freedom which properly describe 
the gravitational field at the quantum level. In view of the successes of classical general 
relativity (GR) the most natural working hypothesis is that they are given by the gauge 
invariant content of the metric tensor g^j^y. Then, trying to construct a quantum field 
theory based upon g^y, a second question arises: what is the (bare) action describing its 
dynamics? The obvious first option to try is the Einstein-Hilbert action underlying classical 
general relativity, but it is well-known that its perturbative quantization leads to a non- 
renormalizable theory. The next logical possibility is to stay within perturbation theory 
but to use a different action. For example, one could try to add terms quadratic in the 
curvature to the Einstein-Hilbert term; this yields a theory which is indeed perturbatively 
renormalizable but nevertheless has to be discarded because of unitarity problems. 

Thus it seems likely that a satisfactory quantum field theory of the continuum metric, 
if it exists, requires us to abandon both the classical action and perturbation theory. Wein- 
berg's "asymptotic safety" scenario P is an approach in precisely this category^ Here the 
basic idea is to define, non-perturbatively, a theory which is based upon a very special bare 
action, namely one which is infinitesimally close to a fixed point r*[5f^i,] of the Wilsonian 
renormalization group (RG). If it exists, the fixed point allows us to take the limit of an 
infinite ultraviolet (UV) cutoff in a controlled, or "safe" way. 

The theory thus constructed will be referred to as Quantum Einstein Gravity or "QEG" . 
^Approaches which (typically) retain the Einstein-Hilbert form of the action but instead abandon the 
traditional continuum include loop gravity |2| and random triangulations 
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Figure 1: The Type Ilia trajectory realized in Nature and the separatrix. The dashed line is a trajectory 
of the canonical RG flow. (Taken from ref. |2(Jj.') 

Using the technique of the exact RG equations in the continuum |1] - jZj a considerable 
amount of evidence has been collected by now |H]-5H| which indicates that an appropriate 
fixed point is indeed likely to exist on the "theory space" consisting of the diffeomorphism 
invariant functionals of gf^,y. Practical calculations typically involve a truncation of this 
theory space. A conceptually quite different approximation leading to similar conclusions 
is the (exact) quantization of the restricted class of metrics admitting two Killing vectors 

uni- 

Most of the investigations using the effective average action formalism |3] - [7j for 
implementing the RG "coarse graining" employ the so-called Einstein-Hilbert truncation 
which retains only Newton's constant G{k) and the cosmological constant A(A;) as running 
parameters. Here k denotes the variable infrared (IR) cutoff, the momentum scale down 
to which the quantum fluctuations of the metric are integrated out. If one introduces the 
dimensionless couplings g{k) = k'^G{k) and X{k) = A{k)/k'^ the flow equations in this 
truncation consist of an autonomous system of two coupled ordinary differential equations: 



kdkg = Pg{g,\) , kdkX = px{g, \). (1.1) 

The beta-functions Pg and Px were first derived in |H], and the resulting equations were 
studied numerically in ref. ^T]. The flow on the g-X-plane displays two fixed points: a 
Gaussian fixed point (GFP) at the origin, and a non-Gaussian fixed point (NGFP) at 
{g^:, A^,) with g^ and A^, both positive. At least according to what can be said on the 
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basis of the Einstein-Hilbert truncation, this NGFP has all the properties necessary for 
asymptotic safety. 

Ref. mi contains a complete classification of all types of RG trajectories occurring 
in the Einstein-Hilbert truncation. Particularly important are those of Type la, Ila, and 
Ilia which, when k is lowered, run towards negative, vanishing, and positive values of the 
cosmological constant, respectively. In ref. [20| the very special trajectory which is realized 
in Nature has been identified and its parameters were determined. This trajectory is of 
Type Ilia; it is sketched schematically in fig. ^ 

For /c — > oo it starts infinitesimally close to the NGFP. Then, lowering k, the trajectory 
spirals about the NGFP and approaches the "separatrix" , the distinguished trajectory 
which ends at the GFP. It runs almost parallel to the separatrix for a very long "RG 
time"; only in the "very last moment" before reaching the GFP, at the turning point T, 
it gets driven away towards larger values of A. In fig. Qthe points Pi and P2 symbolize 
the beginning and the end of the regime in which classical general relativity is valid ( "GR 
regime"). This section of the trajectory is virtually identical to a canonical one for which 
G and A are /c- independent. The classical regime starts soon after the turning point T 
which is passed at the scale /ct ~ 10"^°mpi corresponding to the (macroscopic!) length 
k^^ ^ 10^^ cm.^ 

In ref. [201 it was argued that there starts a regime of strong IR renormalization effects 
to the right of the point P2 which might become visible at astrophysical and cosmological 
length scales. In fact, within the Einstein-Hilbert approximation, trajectories of Type Ilia 
cannot be continued to the extreme IR {k 0). They terminate at a non-zero value 
of k as soon as the trajectory reaches A = 1/2. (Close to the question mark in fig. H) 
Before it starts becoming invalid and has to be replaced by a more precise treatment, 
the Einstein-Hilbert approximation suggests that G will increase, while A decreases, as 

A / 1/2 m]. 

The Type Ilia trajectory of QEG which Nature has selected is highly special or "un- 
^ Since it is only the cosmological constant which shows a significant running between k = mpi and k — 
/ctj and since A cannot probably be measured in a millimeter-size laboratory experiment, one nevertheless 
would not expect to see violations of classical general relativity immediately above kr PH] . 



3 



natural" in the following sense. It is fine-tuned in such a way that it gets extremely close 
to the GFP before "turning left" . In [20] it was shown that the coordinates and At of 
the turning point are both very small: qt = Xt ^ 10^^°. In the GR regime, g decreases 
from g{k) = 10~™ at a typical terrestrial length scale of = 1 meter to g{k) = 10~^^ at 
the solar system scale of = 1 astronomical unit. Extrapolating to cosmological scales 
one finally has g{k) = 10^^^° when k equals the present Hubble constant Hq. 

In this analysis the two free parameters which uniquely characterize a Type Ilia tra- 
jectory where derived from the measured values of G{k) for laboratory values of k, and 
A(/c) at ~ Hq. It can be argued that the IR renormalization effects, if they exist, could 
not change G and A by many orders of magnitude between P2 (at solar system scales, say) 
and cosmological scales. 

Indeed, the present Hubble parameter k = Hq is approximately the scale where the 
Einstein-Hilbert trajectory becomes unreliable. The observations indicate that today 
the cosmological constant is of the order Hq. Interpreting this value as the running 
A{k) at the scale k = Hq, the dimensionless X{k), at this scale, is of the order unity: 
X{Hq) = K{Hq)/HI = C(l). So it is quite precisely near the present Hubble scale where 
the Einstein-Hilbert truncation becomes insufficient for a description of the trajectory Na- 
ture has chosen. 

The "unnaturalness" of Nature's gravitational RG trajectory has an important conse- 
quence. Because it gets so extremely close to the GFP it spends a very long RG time in 
its vicinity because the /3-functions are small there. As a result, the termination of the 
trajectory at A = 1/2 is extremely delayed, by 60 orders of magnitude, compared to a 
generic trajectory where this happens for k near the Planck mass. In ref. j2^ it was argued 
that this non-generic feature of the trajectory is a necessary condition for a long classical 
regime with G,A = const to emerge, and any form of classical physics to be applicable. 

It was also shown that for any trajectory which actually does admit a long classical 
regime the cosmological constant is automatically small. Nevertheless, the fine-tuning 
behind the "unnatural" trajectory Nature has selected is of a much more general kind 
than the traditional cosmological constant problem 01 121|: the primary issue is 

the emergence of a classical space-time; once this is achieved, the extreme smallness of the 
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observed A (compared to comes for free. (See [20] for a detailed discussion.) 

Knowing at least the qualitative features of the RG running of the gravitational parame- 
ters one can try to use this information for investigating how quantum gravity effects modify 
the classical Friedmann-Robertson- Walker (FRW) cosmology. An immediate consequence 
of an approach of this kind is that the cosmological constant becomes a time-dependent 
quantity [22] so that in principle it should be possible to understand its (tiny) value today 
as the result of a dynamical evolution process. In this respect the RG approach has certain 
features in common with the quintessence models [2nil2Z|- The main difference is that the 
evolution is driven by the vacuum fluctuations of the gravitational field itself, and no extra 
quintessence field needs to be introduced. 

In ref. [2H] a first application of the RG flow from QEG to the cosmology of the early 
universe has been described. The tool used there was a kind of gravitational "RG improve- 
ment" [2ni - jinj- The running couplings G{k) and A(A;) where converted to functions of 
the cosmological time t by identifying the RG scale k with 1/t or, what was the same in 
this context, the Hubble parameter H{t). Then the resulting G{t) and A{t) were inserted 
into the Einstein equations for a homogeneous and isotropic universe. In the fixed point 
regime where 

Gik)=gjk\ A{k) = X,k\ (1.2) 

and G{t) oc A{t) oc it was then possible to find exact analytic solutions for the 
Robert son- Walker scale factor a{t). Provided this improvement of the cosmological evolu- 
tion equations does indeed encapsulate the leading quantum effects, the fixed point solution 
should describe the universe at times much earlier than the Planck time since ()1.2|) is valid 
only as long as k ^ mpi, the Planck mass being the lower boundary of the asymptotic 
scaling region. A strong argument in favor of the validity of the approach is that gravity 
becomes weakly coupled for t — > since G{k), vanishing for k oo, is an asymptotically 
free coupling. 

In particular for the equation of state p = p/3 and flat 3-sections {K = 0) the flxed 

point cosmology has various remarkable properties.^ It is completely scale-free; the scale 
•^In ref. l28l also fixed point solutions for the more general equation of state p = wp and (pseudo) spherical 
3-spaces {K — ±1) were obtained. In a different context similar cosmologies were investigated in |24J. Fixed 
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factor and the matter energy density behave as a(t) oc t and p{t) oc respectively. For 
every RG trajectory interpolating between the fixed point behavior ()1.2|1 and the classical 
G,A = const, the fixed point solution is a universal attractor in the space of improved FRW 
cosmologies in the sense that they all approach the fixed point solution for t — ^ 0. Stated 
the other way around: the fixed point cosmology prepares the initial conditions for the 
subsequent cosmological evolution which is essentially classical. Those initial conditions 
are such that the energy density today equals precisely the critical one, pcrit- No fine- 
tuning is necessary here, only a discrete choice, picking fiat rather than (pseudo) spherical 
3-spaces Furthermore, as a result of the linear expansion a{t) oc t at early times, the 
RG improved cosmology has no particle horizon. This is another feature which could be 
of phenomenological importance. 

Up to now the solution to the coupled system of RG- and cosmological evolution equa- 
tions is known only in the fixed point regime. The purpose of the present paper is to 
analyze the complete solution all the way from the Big Bang to scales of the order of 
the present Hubble constant. Using mostly numerical methods we shall discuss in detail 
whether or not the Einstein-Hilbert truncation is sufficient in order to connect the UV 
fixed point regime (Planck era) to the present era of the universe. 

The remaining sections of this paper are organized as follows. In Section 2 we briefiy 
review some properties of the coupled RG- and cosmological evolution equation, point out 
a conceptual difficulty related to the fact that this system is overdetermined for any "rigid" 
cutoff identification, and propose a solution to this problem. Then, in Section 3, we present 
the RG equations to be used in the following. Section 4 contains an investigation of all 
those properties of the RG improved cosmology which do not depend on the detailed form 
of the RG trajectory but rather refiect certain general properties of the "theory space" 
on which the RG fiow takes place. In Section 5 we explore the actual cosmological time 
evolution; in doing so we advocate the point of view that this evolution is primarily an 
evolution with respect to the scale k. Section 6 contains a brief summary and discussion 
of the results. 

point cosmologies of the late universe are discussed in [M\ 02] . 
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2 RG improved Einstein's equation 

Following [2H1 inij we investigate homogeneous and isotropic cosmologies described by a 
standard Robertson- Walker metric containing the scale factor a{t) and the parameter K = 
0, ±1 which distinguishes the three possible types of maximally symmetric 3-spaces of 
constant cosmological time t. The dynamics is governed by Einstein's equation G^^ = 
—Ag^^ + SnGTi^ii, where T^'^ = diag[—p,p,p,p] is a conserved energy momentum tensor for 
which the equation of state p{t) = wp{t) with an arbitrary constant w > — 1 is assumed. 

To start with, we assume that we have fixed a certain "cutoff identification" k = k{t). 
Then every RG trajectory k i— {G{k), A{k)) obtained by solving the RG equations gives 
rise to time dependent functions G(t) = G{k = k{t)) and A(t) = A{k = k{t)). Replacing 
the constants G and A in Einstein's equation with these functions, and specializing for a 
Robert son- Walker metric, the field equations boil down to 

aV K 1 ^ 8n ^ , , 

- + = A+ Gp, 2.1a 
a / 6 6 

p + 3{l + w)-p = 0, (2.1b) 
a 

A + 8npG = 0. (2.1c) 

Here the "dot" denotes the derivative with respect to the cosmological time t. Eq. ()2.1aj) is 
the familiar Friedmann equation, albeit with a time dependent G and A, and eq. (j2.1bjl is 
the conservation law D^T^^ = 0. Eq. ()2.1cp represents an additional consistency condition 
which is necessary for the integrability of Einstein's equation. If ()2.1c|l holds true, its RHS 
has a vanishing covariant divergence, as has its LHS by virtue of Bianchi's identity. 

For G{t) and A{t) fixed, the equations ()2.1al b.c) constitute a system of 3 equations 
for only two unknowns, a{t) and p{t). In general this system is overdetermined and, for 
generic G{t) and A(t), will not admit any solution for a{t) and p(t). This can be seen as 
follows. The time dependence of p follows directly from ()2.1c|l . 

P(t)---^, (2.2) 
' 8n Git) ^ ' 

and integrating ()2.1b|l yields pa?^^'^ = A4/{8tt), with M. a constant of integration. Com- 
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billing this relation with ()2.2|) we find the time dependence of the scale factor: 

l/(3+3«)) 



a{t) 



(2.3) 



A(t)_ 

Already at this point we know a{t) and p(t), but we have not used the improved Fried- 



mann equation ()2.1a|l yet. Checking whether ()2.2|) and ()2.3|) also solve ()2.1a|l . one usually 
discovers that this is not the case: the three equations ()2.1a| b.c) are not integrable for 
generic G{t) and A(t). 

At this point there are two logical possibilities: either the improvement of the field 
equations is not an adequate way of exploiting the RG information, or one can ultimately 
arrive at a description which is self-consistent within this framework by modifying G{t) 
and A(t) in such a way that all three equations ()2.1al b.c) are satisfied. As for the second 
possibility, there are two obvious ways of modifying the time dependence of G and A. On 
can change the cutoff scheme [HI EH], i-e., the operator IZki^—D"^) which suppresses the 
"slow" modes in the path integral, thus slightly changing the trajectory k ^ {G{k), A{k)), 
or one can change the "cutoff identification" k = k{t). The details of both the trajectory^ 
and the cutoff identification are non- universal, i.e., unphysical, but in combining G{k) with 
k = k{t) to obtain G{t) = G{k = k{t)) the non-universalities can cancel to some extent,^ 
provided k{t) is chosen appropriately. 

The length scale set by the IR cutoff^, I ~ /c"^, can be visualized as the variable 
resolution of the "microscope" with which space-time is observed |^ since the effective 
average action describes the dynamics of fields averaged over volumes of linear extension 
£ In cosmological applications it is sensible to readjust the resolution of the microscope 
as the universe becomes larger. In which way precisely I is increased as the scale factor 
grows depends on the interpretation one wants to give to the coarse grained picture of the 
expanding universe. In |^ the ansatz k (xl/t was used, motivated by the fact that when 
the age of the universe is t, no process with frequency smaller than 1/t can have occurred 
yet. Also the Hubble scale k oc H{t) would be a natural choice since in cosmology the 
"^The critical exponents of the fixed points, for instance, are universal (i.e., 7?.fc-independent) features 

of the RG flow, but not the precise shape of the trajectories. 

^See ref. for a particularly transparent example in a black hole context. 
^See ref. ^3] for a more detailed discussion. 
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Hubble length = 1/H{t) measures the size of the "Einstein elevator" outside which 
curvature effects become appreciable. 

In the present paper we shall adopt the following strategy for choosing k{t). In order to 
circumvent the problem of the overdetermined equations ()2.1al b.c) we shall not a priori fix 
the cutoff identification k = k{t) in a. rigid way, but rather derive it from the field equations 
themselves, imposing the condition that those equations should be integrable. On the RG 
side, we shall stick to a fixed cutoff scheme. Then, choosing initial conditions for the 
trajectory, we obtain uniquely defined functions G{k) and A(A;). Next it is important to 
note that, given the fc-dependence of G and A, we immediately know a and p as functions 
oft. 



a{t{k)) 



-M 



G'jk) 
A'{k) 



l/(3+3«)) 



(2.4a) 



On the LHS of these relations t{k) is the functional inverse of the cutoff identification 
k = k{t). The eqs. (^Z^h) follow from (123) and (Q, respectively, since G = (dk/dt) G', 
A = (dk/dt) A', and the derivative of k{t) drops out from the ratio G/A. (Here and in the 
following the prime denotes a derivative with respect to k.) 

With ()2.4a^ b) the second and the third of the eqs. ()2.1a| b.c) are satisfied. We can now 
try to make the whole system consistent by allowing the Friedmann equation ()2.1a|) to fix 
the relationship between k and t. Inserting ()2.4at b) into eq. ()2.1a|) the latter yields the 
following differential equation^ for k = k{t) 



dk a 
dt a' 



A Stt ^ 



(2.5) 



which can be solved (numerically) for any given trajectory.® The RHS of ()2.5|1 is completely 

specified in terms of RG data: G{k) and A(A;) obtained directly from the trajectory, and a 

and p are given by ()2.4a| b). 

This method provides us with a solution {a{t), p{t), k{t)} for any trajectory (^^(A;), A(A;)), 

but it is not clear a priori whether this solution is physically sensible. It can have a mean- 

^Here we have chosen the sign of the square root such that t increases with decreasing k, i.e., late 

cosmological times wih be associated with smah /c-values on the RG-trajectory. 
similar strategy was used in ref. |43| . 
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ingful interpretation only if the continuous readjustment of the "microscope's" resolution 
described by = k{t) is correlated with the expansion of the universe in a transparent 
and, in particular, monotonic way. This will indeed turn out to be the case: We shall see 
that the function k{t) obtained dynamically is reasonably close to k{t) oc H{t) during all 
epochs of the history of the universe. 



3 The RG equations 

The RG equations of QEG in the Einstein-Hilbert truncation are given by (jl.ll) . The 
pertinent theory space is the (7-A-plane coordinatized by the dimensionless couplings g and 
A. The (dimensionless, /c-independent) beta-functions Pg and (3x are the components of a 
vector field on the theory space. In ref. [Hj they were obtained for d space-time dimensions. 

In d dimensions, the dimensionful couplings are given by G{k) = k'^~'^g{k) and A(A;) = 
k'^X{k). Their fc-dependence is governed by the system of equations 



k ^ G{k) = Pg{G, A,k), k^ A{k) = /3a(G, a, k) 
The (dimensionful, fc-dependent) beta-functions Pg and /3a read [Hj: 



PG{G,A,k) =r]NG, 
/3a (G, A,k)=mA + ^(47r)i"'^/2 ^ 



2d{d+l)^'^/,{-2A/k 
- 8rf<l>;^/2(0) - d{d + l)r/^|.i/2(-2 A/k') 
Here tjn denotes the anomalous dimension of y^-R; it can be expressed as 



Vn 



l-gB,{\) 



with the abbreviations 



Si(A) = -{Any-"/' 



52(A) = --(47r 




.l-d/2 



did + l)$^/2_i(-2A) - 6did - l)$^/2(-2A) 
-4rf$i/2_i(0)-24$^_2(0) 
did + l)<l^/2_i(-2A) - 6did - ml/,i-2\) 
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(3.1) 



(3.2a) 
(3.2b) 



(3.3) 



(3.4) 



The "threshold functions" $^ and $^ are given by certain integrals which depend on the 
cutoff scheme, i.e., on RkijP'). In Appendix Elwe discuss them for the cutoff schemes used 
in the present paper: the sharp cutoff jTTj, the exponential cutoff |H| and Litim's optimized 
cutoff jlS]. Note that the system of differential equations (jl.ll) is an autonomous one, while 
dSH) is not. 

In this paper we shall solve (jH.lj) for d = 4 numerically and use the resulting RG 
trajectories {G{k), A{k)) for the improvement procedure outlined in Sectional Since the 
cosmological equations contain matter we should, in principle, use the beta-functions of 
gravity coupled to the corresponding matter system. Since we are interested in a qualitative 
understanding only, we shall use the flow equations of pure gravity, however. As in ref . [20] , 
our analysis is based upon the explicit assumption that the matter fields do not change the 
gross qualitative features of the pure gravity RG flow. Unless we know with certainty what 
the matter fields in Nature are we cannot anyhow decide on a theoretical basis whether or 
not this assumption is really correct. 

4 Cosmology on theory space 

The cosmological evolution with respect to the Robert son- Walker time t is related to a 
scale-evolution via the cutoff identification k = kit). Therefore we may think of the 
history of the universe as a curve in the truncated theory space, k {g{k), A(A;)). 

Remarkably, certain properties of the universe at a given t depend only on the point 
[g, A) in theory space where the universe happens to "sit" at that time, but not on k or 
on the form of the trajectory. Examples include the matter and vacuum energy densities 
(divided by the critical density) and the deceleration parameter. 

The present section is devoted to those properties of the RG improved cosmology which 
are directly related to the theory space (the ^f-A-plane) and can be analyzed without first 
solving for the RG flow. 

There exist two curves on the g-X-p\a.ne which are important for a qualitative under- 
standing of the improved cosmologies: the "fi-line" at which /5a = 0, and a line on which 
1]]^ diverges. 



11 



4.1 The Q-line 

By combining eqs. (j2.4bj) and (jH.lj) we can express the matter energy density directly in 
terms of beta-functions: 

1 ,a^{G{k),A(k),k) 

Let us apply this formula to the class of cosmologies defined by the following properties: 
(i) They are not re-contracting, i.e., at least for t —>■ oo their scale factor a{t) increases 
monotonically with t, and p oc l/a^^^"^ decreases correspondingly, (ii) The beta-function 
Pg does not diverge during the evolution. We shall argue later on that the cosmology 
realized in Nature is likely to belong to this class. 

For these cosmologies we have a{t —>■ oo) cxd and p{t —>■ oo) 0. Hence their formal 
endpoints have p = at "t = oo" due to an eternal dilution of matter. In view of eq. ()4.1|1 
this means that /3\ vanishes for t ^ oo since, by assumption, jSc is always finite: 

limpJGikit)),Aikit)),kit)) =0. (4.2) 

In order to analyze the implication of this condition we first consider the more general 
equation 

/?a(g(A;),A(A;),A;) =0 (4.3) 

where no cutoff identification is invoked and k is considered the independent variable. 
Remarkably, when re-expressed in terms of the dimensionless couplings g = k'^G and 
A = A/A;^, eq. ()4.3|) becomes fc-independent. (Here and in the following we specialize to 
d = 4.) Using eq. ()3.2b|) . we see that /5a = is equivalent to 



r]N{g,\)X + 

ZTT 



10 <1'2(-2A) - 8$^(0) - 5 A) 



. (4.4) 



This provides a condition on g and A. It defines a curve on the (7-A-plane which we shall 
refer to as the "f2-line". It is the locus of all possible endpoints ( "i7-points" ) for the RG 
trajectories k ^-^ {g{k), X{k)) belonging to the cosmologies with eternal expansion and 
Pg < 00. 

What will turn out crucial for the cosmologies of this class is that the fi-line coincides 
by no means with the boundary of the theory space, but is a well defined, albeit scheme 
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Figure 2: The il-line for the optimized cutoff (dashed line), the exponential cutoff with s = 2 (dotted 
line) and the sharp cutoff (solid line). The fi-line prevents cosmological solutions to run into the line 



dependent curve which hes at least partly in its interior. In fig. |21it is plotted for the three 
cutoffs discussed in Appendix A. 

Let us assume we apply the method outlined in Section |21 to a RG trajectory of Type 
Ilia such as the one depicted in fig. ^ Then it can happen that even after an infinite 
cosmological time the RG trajectory on the (7-A-plane does not reach the regime of possibly 
strong IR effects (A < 1/2) but rather stops at some point on the fi-line. In this case the 
f2-line is "screening" the region of strong IR renormalizations. We shall see that, because 
of this mechanism, there are cosmologies which would never experience the conjectured 
strong IR quantum effects and remain essentially classical for t — 00. 

However, the situation is somewhat involved since for most cutoff schemes the actual 
boundary of theory space is not precisely the straight line at A = 1/2, but rather a com- 
plicated curve slightly left of it on which 1]^ diverges already. (The sharp cutoff is an 
exception in the sense that t^at diverges at A = 1/2 only.) Therefore, in order to find out 
whether there can be a cosmological era with non-trivial IR effects, we must know the 



A = 1/2. 
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Figure 3; The lines of diverging rjN for the exponential cutoff with s = 2 (left diagram) and the optimized 
cutoff (right diagram). The solid lines correspond to the locus where rjN diverges, while the dashed lines 
are the corresponding f2-lines. The location of the NGFP is indicated by the dot. 

relative position of the f2-line and the boundary hne with |?7Ar| = oo. This will be the topic 
of the next subsection. 

4.2 Diverging anomalous dimension 

Eqs. ()3.2a| b) show that G and A are strongly renormalized when the anomalous dimension 
7]N is large. This happens close to certain curves on theory space on which rj^ diverges. 
As long as |?77v| is only moderately large we tend to believe the predictions of the Einstein- 
Hilbert truncation. 

The anomalous dimension ()3.3|1 can diverge in two ways: 

(i) The numerator of the RHS of ()3.3|) can diverge, gBi{X) oo. This happens only 
along the line A = 1/2. 

(ii) The denominator on the RHS of ()3.3|1 can vanish, 1 — gB2{X) = 0. 

Analyzing the second possibility for the cutoff schemes from Appendix El we find that, for 
the sharp cutoff, this condition is never fulfilled for g > 0. For the exponential and the 
optimized cutoff, however, it is satisfied along a curve on the ^f-A-plane which lies to the 
left of A = 1/2 (for g > 0). This line can be found numerically; it is displayed in fig. El 

As a result the boundary of the theory space is the straight line A = 1/2 for the 
sharp cutoff, and the curved lines of fig. 01 for the other two cutoffs. On these boundaries, 
r^iv — oo, and no RG trajectory can be integrated beyond them. (See also 46j.) 
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In fig. 01 the solid lines correspond to the locus where 77 at diverges, while the dashed lines 
indicate the position of the fi-line. For the Type Ilia trajectories coming from the left, 
the f2-line generically screens the singularity at A = 1/2, but it does not provide a general 
screening mechanism for the other divergences which arise from a vanishing denominator in 
r/TV. We observe that for a certain value g = g^^t the fi-line and the locus of r/Tv-singularities 
intersect. If g > gint the r^jy-singularities occur at smaller values of A than the f2-line. In 
this regime we find no screening of the ?77v-singularity. 

For g < (yfint, however, the situation reverses and the f2-line occurs at smaller values 
of A than the rjN singularity. In this case, both types of divergences in rj^, the ones 
coming from A = 1/2 and the vanishing of the denominator in 1]^ are shielded for all 
three cutoff schemes. This implies that a Type Ilia trajectory which is sufficiently close 
to the g = 0-axis in the IR will always be stopped at the f2-line before the anomalous 
dimension t]n diverges, independently of the cutoff scheme employed. In view of the tiny 
value g = (9(10~^^°) for A = (9(1/2) this should in particular be the case for the trajectory 
realized in Nature. 



4.3 Energy densities and deceleration parameter 

In order to study the evolution of the energy densities associated to the matter and the 
cosmological constant, it is useful to introduce the relative matter and vacuum energy 
densities: 

We assume a spatially flat universe {K = 0) and therefore have Qtot = + = 1 as 
a consequence of the system ()2.1a| b.c) [28j. Using the Friedmann equation ()2.1aj) . we can 
eliminate the Hubble parameter from the equations for Qm and Q^'- 

1 G'/G 1 A'/A , , 

= 1 — — = 1 — f4 6) 

A'/A ' G'jG ■ ^ ^ 

This motivates deflning the function 

m^^. (4.7) 
15 



which completely determines the evolution of and Q^: 

Y{k) „ 1 



n 



M 



1-Y{k)' " 1-Y{k) 

Here Y{k) is understood as being evaluated along a particular solution of the RG equations. 
The above relations imply that, for Y{k) = 0, the total energy density is completely 
dominated by the cosmological constant {Qm = 0, I^a = 1), while Y{k) = ±oo corresponds 
to complete matter domination {Qm = I^^a = 0). 

These relations can be used to express the deceleration parameter 



a a 



a? 



(4.9) 



in terms of Y{k). From the improved system (j2.1a[ b.c) one derives [S] that, as in the 
classical case, 

q = \{'iw + l)VlM-^K. (4.10) 



Substituting ()4.8jl we obtain 



Y 



3m; + 1 



Y + l 



Writing 



Y{k)=y(g{k),\{k)^ 
in terms of the dimensionless couplings A it turns out that 



y{9A) = 1 + 



9 



10 <1'2(-2A) - 8<1>^(0) - 5 riN^{~2 A) 



(4.11) 



(4.12) 



(4.13) 



2717] N{g A) A 

has no explicit /c-dependence. By definition, F is a function of /c, while ?/ is a function of 
g and A. Hence VLm-, ^a, and q are completely determined by the values of g and A: 

y{gA) 1 l + {^w + l)y{g,X)/2 



M 



q 



l-l/(^,A)' l-y{g,X)' y{g,X)-l 

Thus it is important to know the properties of the function y{g,X). 



(4.14) 



We start by comparing eq. ()4.13|) to the condition determining the f2-line, eq. ()4.4j) . 
This shows that y{g, A) = at the fi-line, so that the energy density of a universe on this 
line is completely dominated by the cosmological constant, = 1, = 0. 
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Figure 4: The function yo{X). Its zero corresponds to the fi-hne at g = 0. The right figure iUustrates 
that near this zero yo{X) is an approximately Unear function. 

In our further analysis we first focus on the g = 0-hmit of y{g, A). This is motivated 
by the fact that the trajectory reahzed in Nature has (? ^ 1 so that this hmit provides a 
good approximation. Setting g = in ()4.13|) we then obtain: 

1 



yo{X) = Urn y{g,X) 



[I0<l>i(-2A)-8<I>^(0)] . 



(4.15) 



27rA5i(A) 

Using the sharp cutoff, this function is displayed in fig. HJ Looking at the left diagram, 
we observe that the domain of definition of y^, (— oo < A < 1/2), contains the following 
special points and intervals: 



A — oo 
-oo < A < 
A = 0+ 
< A < A^_ijng 

^n-line < A < 1/2 
A = 1/2 



-oo 



2/o(A)\l, 

yo(A) > 1 , 

2/o(A) <0, 
< yo{\) < 1 

yoW = 1 . 



(4.16) 



Fig. El shows the full function y{g,X) on the g-X-plane. This figure illustrates that the 
behavior described above also extends to the region where g > so that the discussion 
given below will also apply to a general trajectory with g > 0. In particular, close to the 
f2-line, y is an approximately linear function of A. 

In terms of Qm and Q^, the features of y found above have the following interpretation. 
The region — oo < A < is associated to the IR behavior of RG-trajectories of Type la. A 
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Figure 5: The function y{g,X). The sohd black Une represents the fi-Une where y{g,X) = 0. Note that 
y(5j \ 0) ^ which corresponds to a completely matter dominated universe. 

cosmology evolving along such a trajectory will asymptote to Qa —oo, Qm —>■ +oo with 
the sum + Qm = 1- 

Next we have the separatrix which, in the IR, ends at the GFP g^, = 0, A* = 0. At 
this point i/q^X = 0+) = — oo, so that the cosmology asymptotes to Qm = I^^a = 0. The 
resulting universe will then be completely matter dominated at late times. 

In the region < A < A^_]jj^g, we find the RG-trajectories of Type Ilia, discussed 
in the previous section. In this case the cosmological evolution ends at the f2-line where 
y{g, A) = 0. Hence these cosmologies lead to a domination of the cosmological constant at 
late times: = 0, = 1. 

An interesting property of this latter region is that it splits into two rather different 
parts. In the first part, where y{g, A) increases rapidly, we observe matter domination while 
in the linear regime Q/y dominates. 

In the region Aj^_jjj^g < A < 1/2 we have i^A > 1 and Qm < 0, which is (probably) 
unphysical. 
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For later use we note the value of the function y at the NGFP ((7*, A*): 



y{g,,K) = -l (4.17) 

This result is most easily proven by simply inserting A(A;) = A*A;^ and G{k) = g*/k'^ into 
eq. (gill). By virtue of eq. dHHI), this implies that Q*m = ^1 = 1/2 and q* = (3w - l)/4 
at the fixed point. 



5 Evolution in time and scale 

The proper way of thinking about an RG improved cosmological history is to visualize it 
as a curve on theory space, the RG trajectory, which, by means of eqs. ()2.1al b.c). induces 
both a Robert son- Walker scale factor a{t) and a cutoff identification k = k{t). In this 
section we "switch on" this k- or t-evolution and study, mostly by numerical methods, the 
coupled system ()2.1al b.c) after having first obtained a trajectory {g{k), X{k)), or rather 
{G{k), A{k)), by solving the flow equations (jHH}- Because of their special relevance we 
shall consider trajectories of Type Ilia only. Since the crucial qualitative features such as 
the existence of an f2-line are the same in all cutoff schemes, we shall use the technically 
convenient sharp cutoff throughout. 

5.1 Initial conditions and RG trajectories 

To investigate the properties of the solutions {g{k), X{k)) numerically we choose our initial 
conditions for the dimensionless couplings g, A along the line connecting the turning points 
{gr, At) of the trajectories. (See fig. El) For the sharp cutoff and close to the GFP this 
line is given by [20] 

^t = ^9t- (5.1) 
For the numerics we supplement this relation with 

^init = 1 ) ^init = , Oinit = 1 ■ (5-2) 

This gives rise to a one-parameter family of cosmological solutions which, loosely speak- 
ing, are characterized by their distance to the GFP at X = 0, g = 0. By definition, the 
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Figure 6: RG-trajectories with the initial conditions given in H5.3|) . The dotted line is the locus of the 
turning points which separate their UV and IR parts. The two almost vertical lines are the q = 0-locus 
for w — 1/3 and the fi-line. 

trajectories pass the turning point at the scale kx, i.e., gi^kx) = Qt and X^kx) = Xt- Thus, 
adopting the relations ()5.2|) amounts to expressing all dimensionful quantities in terms of 
appropriate powers of kx- (In this parameterization the Big Bang occurs at t < 0.) 
Starting with these initial values we evolve the solutions in two directions: 

• Towards the UV {k > kx) where the RG-Flow will be attracted towards the NGFP. 

• Towards the IR {k < kx) where the cosmology will run towards the f2-line. 
For explicitness we consider the three sample solutions with 



These examples are sufficient to illustrate the general trend for decreasing qt and to un- 
derstand the qualitative properties of the RG-trajectory realized in Nature which has an 
extremely tiny gf- The RG-trajectories obtained by numerically solving ()3.1|) with these 
initial conditions are shown in fig. IHl We see that decreasing results in squeezing the 
solution into the corner between the separatrix (connecting the NGFP and the GFP) and 
the g = 0-axis. In fact, within the resolution of fig. IHlthe two trajectories with gx = 10^^ 
and gT = 10~^ run virtually on top of the separatrix (for k > kx) and the horizontal axis 
(for k < kr)- 



gT = io 



-1 



gx = 10 ^ , and gx = 10 



-3 



(5.3) 
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All three trajectories are of Type Ilia and terminate at the boundary A = 1/2. For 
gr = 10^^, 10~^, and lO^'^ the termination scales are found to be fcterm/^T = 0.178, 0.0612, 
and 0.0196, respectively. These numbers are consistent with the rough estimate fcterm/^T = 
Oi^^Jgr) derived in [201 They confirm that the closer the trajectory approaches the GFP 
the later (in fc) it terminates. 

5.2 UV vs. IR branch 

The two branches of the trajectory, k > kT and k < kx, give rise to two corresponding 
branches of the cosmological solution which we shall refer to as the "UV cosmology" and 
"IR cosmology", respectively. In the UV cosmology we are going to use the equation of 
state with w = 1/3 ("radiation"), while we employ w = ("dust") for the IR cosmology. 
In reality the universe has passed the scale k = kx well inside the radiation dominated era, 
but since we are mostly interested in qualitative features of RG improved cosmologies this 
aspect is inessential. 

Next we shall apply the method described in Section |21 to the trajectories of the previ- 
ous subsection and discuss their UV- and IR-branches in turn. In all examples the Hubble 
parameter is calculated using eq. ()B.7|) derived in Appendix |Bj and the deceleration pa- 
rameter is determined with the help of eq. ()4.14|1 . 

5.3 RG-improved UV cosmologies 

The numerical results for the UV cosmologies resulting from the (^r-values ()5.Hj) are dis- 
played in figs, m - [ini Fig. [7| shows the scale factor a(t), and fig. |H1 the corresponding 
Hubble constant H = a/a and deceleration parameter q. In fig. El Newton's constant 
G{t) = G{k{t)) and likewise the cosmological constant is plotted as a function of time. 
Here t is measured in units of k^^, and the origin t = corresponds to k = kx- 

^The proportionality constant is fctcrm/^T ~ 0.612y/gT- The first case gx — 10~^ does not quite fit into 
the pattern, since this trajectory is not yet sufficiently deep in the GR regime. 
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Figure 8: H{t) and q{t) for the UV cosmologies starting with = 10 ^ (short dashes), gr = 10 ^ (long 
dashes) and gr — 10~^ (sohd Une). 

5.3.1 The regimes of the UV cosmologies 

We observe that the cosmologies possess an initial singularity ("Big Bang") at a time 
ts < where a{tB) = and H(t \ ts) — > oo. For t > we can distinguish the following 
regimes: 

(A) The NGFP regime 

Immediately after the Big Bang, for k 3> mpi, the RG trajectory is well approximated by 
the constant functions g{k) = and X{k) = A*. They correspond to the fc-dependence 
fll.2|) for which the system ()2.1al b.c) can be solved analytically. This yields the fixed point 
solution found in ref. [2Hj. For w = 1/3 it yields a oc (t — t^), H = l/(t — ts), q = 0, 
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Figure 9: The functions G{t),A{t) for the UV cosmologies arising from the initial conditions = 10 ^ 
(short dashes), gr — 10^^ (long dashes) and gx — 10~^ (solid line). 

G oc (t — tfi)^, and A oc l/{t — tsY- The numerical solutions do indeed show this behavior 
close to ts- In particular, G vanishes and A diverges at the Big Bang. 

(B) The linear regime of the NGFP 

At slightly later time, corresponding to smaller scales k, the RG flow can be linearized 
about the NGFP CH]- The linearized trajectories {g{k), X{k)) are spirals about the 
NGFP, characterized by two critical exponents 6' and 6" [1T\ ITU] . The most prominent 
cosmological signature of this regime is the oscillatory behavior of q{t) seen in fig. |Hland, 
in a slightly magnified way, in fig. ^1 Directly at the NGFP one has g = 0; when the 
trajectory starts circling about the NGFP there are phases with both g > (deceleration) 
and g < (acceleration). In fact, there exists an infinite sequence of phases with either 
sign. When one approaches the Big Bang from above {t \ ts), Q decreases from about 
q ^ 1, has a first zero, becomes negative, has another zero, becomes positive for a period 
shorter than the previous one, then again becomes negative for a time span shorter than the 
previous one, and so on. The Big Bang is approached by an infinite sequence of oscillations 
in q, of decreasing amplitude and decreasing duration. 

(C) Crossover and linear regime of the GFP 

At still lower scales, or later times, the RG trajectories leave the linear regime of the NGFP 
and very quickly "cross over" towards the GFP. Linearizing about the GFP, the resulting 
cosmologies are again easily understood analytically. One finds [201 that approximately 
G{k) = const and A(A;) = A(A;'p) [1 + {k/k'rY]/2. If k is related to t in a monotonic way. 
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Figure 10: Oscillatory phases of the deceleration parameter for the UV cosmologies starting with 
qt — 10^^ (short dashes) and = 10^'^ (solid line). 

we expect G to stay constant and A to decrease towards its value at (where it becomes 
constant, too). This is exactly what we observe in fig. 01 

In the linear regime of the GFP, for the time t < 0, i.e., before the turning point of 
the trajectory, the cosmological constant is not important for the three solutions. At least 
for gx small they have ~ and ~ 1 there, yielding g ^ 1. This is exactly the 
plateau-value of q{f) which is approached in the second plot of fig. |H1 

5.3.2 Dynamical determination of A; = A;(t) 

As we explained in Section |2l our method determines the cutoff identification k = k{t) 
dynamically in such a way that the field equations are integrable by construction, the 
risk being that the k{t) thus obtained does not lend itself for a clear interpretation of the 
resulting coarse-grained cosmology. In the first diagram of fig.^Jwe display the relationship 
between k and t by plotting {t{k)k)~^ as a function of \n{k/kj)- The motivation for 
this presentation is as follows. If this function is fc-independent, t{k)k = const, we have 
t{k) oc 1/k or vice versa k{t) oc 1/t which is the rigid cutoff identification used in [2H1E]- 
In order to get the {t{k)k)^^ vs. ln(A;/A;T) plot, we have shifted t{k) such that the Big Bang 
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Figure 11: it{k)k) ^ and H{t{k))/k for the combined UV and IR cosmologies with gx = 10 ^ (short 
dashes), gr = 10^^ (long dashes) and gx = lO^'' (solid line). 

{k = oo) corresponds to t{k = oo) = O}^ In an analogous fashion the second plot in fig. 
im shows H{k)/k as a function of k; this function would be constant for the identification 
k{t) oc H{t). 

From fig. ^2 we see that both k oc H{k) and k oc l/t{k) provide a valid approximation 
to the true cutoff-identification in the UV domain {k > kx) if one is interested in the 
overall qualitative features only. In this domain the approximation k oc l/t{k) performs 
somewhat better than k oc H{k). 



5.3.3 Oscillatory inflation 

Finally we return to the phenomenon of the g(t)-oscillations in the NGFP-scaling regime 
which we described in 15.3.11 (B). In order to analyze the sign-flips of the deceleration 
parameter we recall from ()4.14|) that q = q{g, A) can be regarded as a function on theory 
space, the t-dependence arising by inserting the trajectory with k = k{t), i.e., q(t) = 
q{gikit)),X{k{t))). 

As for the function q{g, A), there are regions on the (7-A-plane where it is positive and 
others where it is negative, the common boundary being the "g = 0-line" shown in flg. IHl 
By virtue of ()4.14j) this line is given by the implicit condition 1 + {3w + l) y{g, A)/2 = 0. To 
obtain the diagram in flg.lHlthe equation for w = 1/3, y{g, A) = —1, was solved numerically. 

^"This is different in figs-IHlandElwhere the UV cosmology corresponds to i < 0. At t <^ fc = fcy = 1, 
the expression [t[k)k)~^ develops a pole, which is avoided by shifting the Big Bang to i = 0. 
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To the right (left) of the q = 0-hne one has g < (g > 0). 

Now it is important to observe that precisely for w = 1/3 the q = 0-line runs exactly 
through the NGFP. In fact, according to ()4.17|1 the NGFP has y{g*, A*) = —1, i.e., it lies 
on the q = 0-line pertaining to if = 1/3. As a consequence, when the RG trajectory spirals 
around the NGFP it crosses the q = 0-line an infinite number of times in either direction. 
This explains the g-oscillations about g = found in the numerical solutions. 

Thus the RG-improvement predicts epochs of accelerated expansion ( "inflation" ) in the 
early universe and one might wonder whether there can be any relation to the traditional 
models of cosmological inflation. To analyze this question we compute the number of "e- 
folds" the universe expands during the periods with g < 0. We consider only the last such 
period because it has the longest duration and the most negative q{t), see fig.HUl Note that 
even there q{t) never becomes as small as —1 which would correspond to a de Sitter- like 
phase. 

We first derive a general formula which can be used to calculate the number of e-folds 
of expansion, N, occurring along any RG-trajectory. This formula will then be used to 
determine for the last period of accelerated expansion occurring before the solution 
leaves the linear regime of the NGFP (see fig. Ej). 

Our starting point is the usual definition of N, 

a{tc 



iV = In 



a{ti 



(5.4) 



where tb and tg denote the cosmological times where the acceleration begins and ends, 
respectively. Using ()2.4aj) we can relate N directly to the RG-trajectory: 

1 



N 



3 + 3w 



In 



(5.5) 



Here and in the following the subscripts b, e denote the corresponding quantity at the begin- 
ning and end of the acceleration period. Eq. ()5.5p can be evaluated further by re-expressing 
the quantities G', A' using eq. ()B.4|) . Writing the result in terms of the dimensionless cou- 
phngs g, X we find 



N 



3 + 3w 



In 



'h' 



+ 



3 + 3w 



In 



A{9h, At 



(5.6) 
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Tp 

' a 








10-1 


2.505 


1.642 


-1.196 


-1.390 


0.901 


10-2 


3.710 


2.848 


-1.196 


-1.392 


0.900 


10-3 


4.886 


4.004 


-1.195 


-1.392 


0.900 


10-1° 


12.811 


11.945 


-1.196 


-1.392 


0.900 



Table 1: Number of e-folds obtained in the last period of accelerated expansion when a Type Ilia 
trajectory is encircling the NGFP. 



where we define 

A{g,X)=g 



^+2;^(-'° ln(l-2A) + 2C(3)-|,-) 



-1 



(5.7) 



In terms of tfie RG-time r = \n{k/kT), and for w = 1/3, this becomes simply 



AT = Tb - Te + - In 



A{ge, Ae) 



(5. 



.^(s-b, Ab) 

Looking at fig. IHl we see that the coordinates of the last and the second to last g = 0- 
crossing, {g^, Ac) and {g^, Xh), are to a very good precision equal for all Type Ilia trajecto- 
ries, so the last term in (j5.8j) gives rise to an almost "universal" contribution. 

We then evaluate N along the last period of accelerated expansion, the semi-circle in 
fig. ini The corresponding values are summarized in table H From this table we observe 
that y4b and Ae are indeed independent of gr as suggested by fig. IHl Furthermore, we find 
that while the values of rt, and increase monotonically when decreasing g^, the difference 
Tb — Tc, too, is independent of the initial condition for gx- This implies that all trajectories 
give rise to essentially the same number of e-folds: ~ 0.9. This number is far smaller 
than the 60 e-folds occurring in the "usual" inflationary scenarios, so that there is clearly 
no direct correspondence. 



5.4 RG-improved IR cosmologies 

In this section we investigate the "IR cosmologies" related to the k < branch of the RG 
trajectories. For the numerical illustration we shall use the initial conditions ()5.3|1 again, 
this time integrating the RG equations towards smaller k. 



27 



. 1 f 
0.05 - 

0.01 - 
0.005 - 



0.001 
0.0005 



0.2 0.4 0.1 



0.8 ^ 1 



0.2 0.4 0.6 0. 



Figure 12: The functions G{k) and A{k) arising from the initial conditions gr = 10 ^ (short dashes), 
gr = 10~^ (long dashes) and gx = 10~^ (solid). 

5.4.1 The running of G and A 

In fig.^Jwe have plotted G{k) and A(A;) for the three sample trajectories. Going downward 
from = 1 (in /cT^-units, as always) the cosmological constant has a weak scale dependence 
at first but then assumes a constant value over a wide range of fc-values. For G, this 
plateau behavior has set in before the turning point even. These plateaus correspond to 
the classical GR regime where, by definition, G,A^ const. At even smaller k, G{k) starts 
increasing until it finally develops a vertical tangent at the termination scale fcterm- (For 
the real RG trajectory with qt ~ 10^^° the increase is less abrupt than it appears in fig. 
1121) For k > kterm, as long as G{k) is not too different from its classical plateau value, the 
Einstein-Hilbert truncation should be reliable still. It is natural to ask, therefore, whether 
the IR increase of G{k) leads to cosmological consequences, as long as we can trust the 
approximation. (Note that A(A;) remains regular at the termination point; A(A;term) = 1/2 

entails A(A;tcrm) = ^term/2-) 

5.4.2 The cosmological history 

Applying the method of Section 2 to the functions G{k) and A(A;) of fig. Elwe obtain the 
cosmological solutions displayed in the figs. Uni (scale factor), HH (Hubble and deceleration 
parameter). 1151 (time dependence of G and A), and 1161 (scale-time relationship). 

The gross features of all three of the typical cosmologies considered are as follows. Near 
k kx the universe is matter dominated {VLm ~ 1, f^A ~ 0) and has a deceleration param- 
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Figure 13: IR cosmology: The scale factor a{t) arising from = 10 ^ (short dashes), (/t = 10 ^ (long 
dashes) and gr = 10~^ (solid line). 
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Figure 14: IR cosmology: H{t) and q{t) for the cosmologies starting with qt — 10~^ (short dashes), 
gr — 10^^ (long dashes) and gx — 10~^ (solid line). Decreasing gx increases the cosmological time-span 
passing until the Hubble parameter becomes constant and delays the transition to a universe dominated 
by the cosmological constant. 

eter q ^ 1/2 therefore. According to fig. El this value decreases towards g = — 1 at asymp- 
totically late times. This asymptotic regime is completely A-dominated: Qm ~ 0,^a ~ 1- 
Decreasing the cosmological time passing during this transition from matter to vacuum 
dominance increases. The value q ~ —0.55, for instance, which roughly corresponds to the 
deceleration parameter observed today, is reached at increasingly late cosmological times 
as qt is becoming smaller. The Hubble parameter approaches a constant value for t —>■ oo, 
indicating that the space-time approaches de Sitter space asymptotically. The kink in the 
a(t)-curves of fig. IT^ marks the onset of the de Sitter behavior. 
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Figure 15: IR cosmology: G{t) and A(t) arising from the initial conditions gx = 10~^ (short dashes), 
gr = 10~^ (long dashes) and gx = 10~^ (solid line). Note the logarithmic scaling of the vertical axis. 
Decreasing gx decreases the asymptotic value of A(t) and G{t) for t ^ oo. 

The transition form matter to A-dominance occurs in classical cosmology as well. To 
what extent are the improved cosmologies modified by quantum gravity effects? Fig. IT5l 
shows that A(t) is decreasing as long as k is not too far below fcy, while G{t) is essentially 
constant. For t > 5A;^^, say, both A(t) and G{t) attain constant values which remain 
unaltered for t ^ oo (at least within the resolution of fig. [T3j) . In the next subsection 
we shall demonstrate that the asymptotic values of A and G are exactly those predicted 
by the fiow equation linearized about the GFP provided qt is small enough. This implies 
that, for qt sufficiently small, i.e., for trajectories getting close to the GFP, there are no 
significant renormalization effects at late times. 

5.4.3 Asymptotic vs. laboratory values and the cosmological constant problem 

Let us come back to the plots of G{t) and A(t) shown in fig. ^1 This figure illustrates 
that decreasing decreases the asymptotic values of A(t) and G{t) for t —>■ oo. Table El 
summarizes these asymptotic values of G and A for various choices of qt- This table also 
contains the corresponding "laboratory" values Giab and Aiab- Their interpretation is as 

follows i2ni- 

Provided the trajectory gets sufficiently close to the GFP it contains a long classical 
regime. (In fig. ^ this "GR-regime" corresponds to the segment of the trajectory between 
the points Pi and P2.) In this regime the dimensionful quantities G and A are approxi- 
mately /c-independent. Linearizing about the GFP one finds for their values, expressed in 
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Figure 16: The function t{k) arising from the initial conditions gr — 10^^ (short dashes), gr = 10^^ 
(long dashes) and gr = 10^'^ (solid). The cosmological time t increases monotonically with decreasing k 
and diverges as the RG trajectory reaches the Sl-line. 
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Table 2: The asymptotic and laboratory values of G and A in units of kr for various choices of gr- 



terms of fcr, 

Aiab = ^At4 = ((P2/47r)^r4 • (5-9b) 

In table El we compare the asymptotic values of G and A to their "lab" values computed 
from ()5.9a| b). We find that for qt of the order 10-2 or smaller the plateau values at 
asymptotically late times virtually coincide with those measured in a "laboratory" in which 
classical GR is known to apply. At qt = 10-^, instead, Gasym is slightly larger than Giab 
while Aasym is slightly smaller than Aiab- 

The interpretation of this result is as follows. For qt sufficiently small, the universe does 
not experience significant IR renormalization effects; in the limit t — > oo it basically keeps 
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Table 3: The asymptotic and laboratory values of G and A in Planck units for various choices of gr- 
Decreasing gx decreases both the asymptotic and the laboratory value of the cosmological constant. 

its values of G and A from the GR regime. Only for comparatively large values of qt we find 
deviations between the GR-formulas ()5.9al b) and the numerically computed asymptotic 
values. Those deviations occur because the RG trajectory does not get sufficiently close to 
the GFP to give rise to a long classical regime where ()5.9a| b) would apply. 

The conventional Planck mass is defined in terms of Newton's constant measured in a 
classical "laboratory": mpi = l/VCiab- This definition, together with ()5.9ajl leads to the 
following important relation between mpi and kr'- 

= Qt iTT'pi ■ (5.10) 

Obviously (7t ^ 1 leads to a large hierarchy kx/mpi ^1. In table El we use ()5.10|) in order 
to express the asymptotic and laboratory values of G and A in terms of the more familiar 
Planckian units. 

As for the cosmological constant problem, the crucial point to be noted is that Aiab ~ 
Aasym is supprcssed by a factor oc relative to its "natural" value mpj. As was discussed 
in detail in ref. fine-tuning the RG-trajectory in such a way that it spends a long 
RG time near the GFP by picking an "unnaturally" small value of qt leads to a long 
(in fc-units) classical regime on the trajectory. Once this is achieved, the solution of the 
cosmological constant problem is for free: 

Alab/"^P1= (<^2/47r)^2 < 1. (5.11) 

The numbers in table El illustrate this important phenomenon. The smaller is Qt, the 
closer the trajectory approaches the GFP, the smaller is the cosmological constant in the 
classical GR laboratory further downstream on the trajectory. 
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The relative importance of the cosmological constant can be summarized as follows. In 
the very early universe, in the NGFP regime, one has Qm = = 1/2 showing that the 
vacuum and matter energy densities drive the expansion on an exactly equal footing. Near 
the GFP, in particular in the classical GR regime, one automatically is lead to ~ and 
Qm ~ 1, a conventional purely matter dominated (flat) FRW cosmology. At infinitely late 
times A takes over again and with Q/y ^ 1 and Qm ~ a de Sitter like behavior sets in. 

5.4.4 fi-line mechanism and dynamical cutoff identification 

In fig. ^Jwe saw that G{k) and A{k) show a strong IR running for small k. It therefore 
comes as a surprise perhaps that G(t) and A(t), according to fig. Cil show no sign of such 
IR renormalizations for t ^ oo; according to the previous subsection, the late cosmology 
is essentially classical. 

This apparent contradiction is resolved in fig.^Jwhich displays the relationship between 
the scale k and the cosmological time t. As expected, t = t{k) is a monotonically decreasing 
function of k. However, an infinitely old universe ("t = oo") does not correspond to A; = 0, 
but rather to a non-zero asymptotic scale 

k^y^ = k{t ^ oo) > . (5.12) 

The function t = t{k) has a singularity at this scale: t{k \ fcasym) oo- Stated differently, 
the inverse function k = k{t) is bounded below by fcasym and, even for arbitrarily late times, 
does not reach arbitrarily small scales: k{t) > /casym for all t. 

The absence of visible IR effects in G{t) = G{k = k{t)) and A{t) = A{k = k{t)) is 
explained by the fact that, for the examples considered, /casym is much larger that the scale 
kterm near which G{k) and A(A;) get renormalized strongly. This is precisely the f2-line 
mechanism we discussed already earlier: for t oo the dimensionless g{k{t)) and X{k{t)) 
approach a point on the f2-line. The endpoint of the cosmology, "a = oo" or p = 0, 
corresponds to this point on the g-X-pla.ne; it is still far away from the \riN\ = cxD-line close 
to which the IR renormalizations would become strong. Therefore, once the universes have 
entered the GR regime at about t ^ S/c^^, say, they remain classical for all later times. 
In this sense the f2-line separates the cosmologically accessible parts of theory space from 
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those with strong IR running. 

The asymptotic de Sitter phase of the universe corresponds to the approximately time 
independent scale k(t) ^ const = /casym ^ herm- 

Let us look more closely at the cutoff identification k = k{t) which was generated 
dynamically by our system of equations. In a slightly different presentation it is displayed 
as the k < /c^-branch in fig. ^2 In this figure {t{k)k)~^ and H(t{k))/k are plotted vs. 
ln{k/kT). If the actual (obviously rather complicated) cutoff identification was close to 
k (X 1/t and k oc H{t), respectively, those plots would show a constant function for all k. 

It is clear, and also confirmed by the plots, that in the IR the real k{t) is quite different 
from k oc 1/t, the reason being that the real k{t) is constant, equal to fcasym, while 1/t 
decreases monotonically for t ^ cxd. 

On the other hand, H{t{k))/k is seen to be approximately /c-independent within a 
factor of less than about 2. Thus, at least at a qualitative or "semi-quantitative" level, 
the actual cutoff identification can be approximated by A; oc H}^ This identification would 
indeed associate a constant scale /Casym to the de Sitter final state, proportional to the time 
independent ifasym = V^asym/3. 

From these observations it is clear how the cosmologies are to be interpreted. Since 
approximately k (x H the length scale characterizing the averaging volume is the Hubble 
radius In = . The cosmological parameters computed correspond to averages over 
the volume The scale ^H characterizes the radius of curvature of the four-dimensional 
space-time, the size of the "Einstein elevator". The asymptotic de Sitter phase makes it 
particularly clear that the temporal proper distance to the Big Bang, t, would not lend itself 
as a coarse graining scale within the improved equations approach. For the UV cosmology, 
in particular near the NGFP [28 , both identifications are equivalent, though. 

5.5 The RG-trajectory realized by Nature 

In ref. [20] the observational (supernova, CMBR, etc.) data were interpreted under the 

assumption that the gravitational RG trajectory which Nature realizes has the qualitative 

features of a Type Ilia trajectory of QEG. It was found that the turning point of this 
^^As an ansatz, k oc H has also been employed in refs. 02 ^ different context. 



34 



trajectory is extremely close to the GFP, and that it is passed at a scale very far below 
the Planck scale: 

gT = O{10~^°), At = 0(10-6°), kr = O{10-^^mpi) . (5.13) 

These order of magnitude estimates do not depend on whether the late universe is, or is 
not affected by IR renormalization effects. The analysis shows that even if there was no 
fi- mechanism the data would not afford a ratio Gcosmoiogicai/Ciab larger than about one 
order of magnitude [20] • 

The universe passes the turning point a fraction of a second after the Big Bang, at a 
Hubble radius of the order in ~ ^ 10^°/pi ~ 10^'^ cm. After that g{k{t)) decreases 
from 10^6° to its present value (7today of about 10~^^°. 

Since g{k) ^ 1 during most epochs of the cosmological evolution, in particular in the 
late (IR) universe, we may use the — * limit of the corresponding beta-functions there. 
This allows for a simple determination of the point on the (^-A-plane where the universe 
resides today. Assuming we know VL^q = f2A(i^today) we can solve the second equation of 
()4.14|1 for the present value of A, Atoday: 

l-2/(0,Atoday)=^^Ad- (5-14) 

Here we approximated (^today ~ 0. Note that the precise value of Atoday is scheme dependent. 

For a rough estimate of Atoday we may assume that the classical FRW cosmology is 
essentially correct at late times. Analyzing the observational data within this framework 
yields Qao ~ 0.7. For this value, and the |/-function corresponding to the sharp cutoff, we 
find 

Atoday ~ 0.320 . (5.15) 

Looking at fig. El confirms that the point ((^today, Atoday) ~ (10^^^°, 0.32) indeed lies to the 
left of the Q-line for the sharp cutoff, as it should be. 

5.5.1 The age of the RG-universe 

As a check of our beta-function formalism, and in order to confirm the conclusions of 
subsection 15. 4| we next derive a relationship between the age of the universe, ttoday, the 
cosmological constant Aiab ~ Aasym, and the measured ^ao- 
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First we use the value of ^ao in order to solve ()5.14p for Atoday(^Ao)- Next we rewrite 
(the reciprocal) eq. ()2.5|1 in terms of the dimensionless couplings g = k'^G, X = A/k"^ and 
expand its RHS up to the leading order in g: 



dt 
dk 



VSk- 



+ Oig) . 



[Bry vrA {BTY - (10 ln(l - 2A) - 4C(3) + ^) 



5f 



ttA fif - (5 ln(l - 2A) - 2C(3)) 



(5.16) 



Here 



3tt 



A 



1^ 



(5.17) 



1-2A (1-2A)2_ 

This expression can be converted to a differential equation for dt/dX by exploiting that in 
the GR regime we have the relation A = Auh/k"^- Its use is legitimate here since the age of 
the universe for trajectories with sufficient "squeezing" is dominated by their GR regime. 
Substituting it into ()5.17|) we obtain the desired general relationship: 

1 



''today 



-^(Atoday) • 



(5.18) 



lab 



Here /(Atoday) is given by the dimensionless integral 



/(A 



today ) 



■^today 



dX 



x/3 



X 



(5f)' 



2A(l-|/o)V2 
ttA (5f )' - (10 ln(l - 2A) - 4C(3) + ^) 



(5.19) 



Bf ttA 5f - (5 ln(l - 2A) - 2C(3)) 

As an example, we evaluate this integral numerically for the Atoday given in eq. ()5.15|) . 
We find 

/(Atoday = 0.32) = 4.19 (5.20) 

Using this value in ()5.19|) . along with y/Ai^ ~ a/ Aasym ~ 10~^°mpi we obtain the age 
^today = 4.19 X 10^° tpi. As cxpectcd, this is essentially the same age one obtains from the 
classical FRW equations if one uses the same input data. 



5.5.2 The deceleration parameter 

We now calculate q{t) along the trajectory realized in Nature below the turning point. As 
in the previous subsection we will work in the approximation g \ 0. We parameterize 
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Figure 17: The deceleration parameter q{t) along the RG-trajectory realized in Nature. The time along 
the horizontal axis is given in units 10^°tpi. At present, t « 4.2, corresponding to q « —0.55. 

the trajectory by its A-values A G [A^, An-une], i.e., k is thought of as a function of A, 
with the inverse A = A(A;). In this parameterization we can use ()5.18|) to obtain t{X{k)) = 
i / VM^I {X{k)) , the cosmological time passing during the evolution along the trajectory, 
while the corresponding deceleration parameter q{X{k)) is obtained from eq. ()4.1H1 . The 
function q{t) along the trajectory can then be obtained as the parametric curve A ^ 



The resulting q{t) is shown in fig. El Note that the time t displayed on the horizontal 
axis is in units lO^'^tpi, indicating that the "crossover" from g(0) = 1/2 to q{t —>■ oo) = — 1 
indeed occurs at cosmological time spans of the order of the age of our present universe. In- 
deed, the present universe corresponds to the age ttoday ~ 4.2 x 10^°tpi and the deceleration 
parameter gtoday ~ —0.55. 

6 Summary and Discussion 

In this paper we presented a coherent renormalization group based framework which al- 
lows for the inclusion of potential quantum gravity effects during the entire cosmological 
history, from the epoch after the initial singularity to the scales of the present universe. 
We demonstrated that the strategy of RG improving the field equations is physically viable 
during all stages of the cosmological evolution. While from a mathematical point of view 



(t(A),g(A)), A G [AT,An-iine]. 
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it is always possible to fix the cutoff scale in a way which renders the modified Einstein 
equations consistent, the result of the dynamical determination of k{t) is not a priori guar- 
anteed to lead to a useful interpretation of the coarse grained cosmology. However, we 
found that the dynamical cutoff identification corresponds to a time dependent averaging 
scale i{t) which, at any time, is approximately equal to the Hubble radius init)- Therefore 
we may conclude that the "microscope" whose observations are described by the RG im- 
proved cosmology has exactly the right "resolving power" to see the large scale structure 
typical of cosmology. 

One of the most interesting features of the early cosmology was the discovery of a kind of 
"oscillatory inflation" . Running time backwards, the initial singularity is approached by an 
infinite sequence of increasingly short time intervals during which the universe accelerates 
and decelerates, respectively. This behavior is due to the fact that the non-Gaussian 
fixed point has complex critical exponents. The oscillatory approach of the Big Bang is 
reminiscent of the classical Belinskii-Khalatnikov-Lifshitz scenario jl?] but contrary to the 
latter the quantum effect found here occurs even in isotropic universes. 

The most interesting aspect of the improved cosmologies at late times is the "fi- 
mechanism" which we investigated in detail. At small scales the RG trajectories of the 
relevant type leave the domain of validity of the Einstein-Hilbert truncation, and it has 
been argued that in the IR there could be a regime where the RG effects become strong 
again. Remarkably, even if this regime exists, the cosmologies we found would never get 
affected by it. Even at arbitrarily late times they remain essentially classical since the 
underlying RG trajectory is never probed beyond a certain "Jl-point" at which rj^ is very 
small still. 

At this point we must ask whether this shielding of this r/Tv-divergence is a reliable 
prediction. Does the universe as a whole really never enter the regime where is large? 
Our argument was based upon two crucial facts: (i) The consistency condition for the 
RG improved Einstein's equation has the simple form ()2.1cj) which leads directly to the 
relationship p oc K'{k)/G'{k) of eq. ()2.4b|) : {ii) The RG flow is such that the zero of ^'{k) 
occurs at a higher value of k than the divergence of G'{k). As a consequence of (ii), p 
vanishes and a oc p^i/(3+3«') caches infinity before |?7Ar| has become large. 
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As for the vahdity of (ii), even though we analyzed the Einstein-Hilbert truncation 
only, this prediction was perfectly stable under a change of the cutoff scheme. Because of 
this robustness, it is likely to be actually correct and not a truncation artifact. 

The much more subtle issue is [i) which is not related to the RG flow per se but to 
the improvement procedure. It is important to note that, at least in this particular form, 
the zero of A'{k) implies an infinite scale factor only in the approach of improving the field 
equations. It would not happen if one improves the solutions of the classical equations by 
replacing G ^ G(t), A i— A{t) there, or if one improves the action functional. The latter 
approach has been investigated in detail in ref. [2Bj- Performing the substitution under 
the space-time integral of the Einstein-Hilbert action one obtains a different consistency 
condition: 



a oo. The improved equation and action approaches are (almost) equivalent in the UV, 
but can lead to different predictions in the IR. There they effectively represent certain 
non-local terms of the ordinary effective action Tk=o generated during the RG running but 
not taken into account explicitly [SHj- The additional terms in which are not present 
in ()2.1c|) are due to the fact that in the improved action approach the scalar fields G and 
A carry energy and momentum which is neglected in the improved equation approach. 

For this reason it cannot be excluded for the time being that the shielding of the 77 at- 
divergence is an artifact of the improvement scheme since, at least via the mechanism 
discussed above, the shielding at A' = occurs only if one neglects the energy and mo- 
mentum carried by the scalars G and A. Clearly more work is needed in order to settle 
this issue. It should be emphasized, though, that whatever is the final answer about RG 
effects in late cosmology it bears no simple relationship to the conjectured modifications of 
gravity at galactic scales. The corresponding cutoff identifications mimicking the essential 
IR-relevant invariants in both cases are extremely difficult to guess beforehand, and in 
particular in a Lorentzian setting where the actual cutoff is (a generalization of) "virtual- 
^^This form of the consistency condition corresponds to the choice 9^i, — 0, see ref. [HHl- 




(6.1) 



We observe that, in this framework, A oc A' 



by no means implies that p 



and 
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ity" rather than "inverse distance" one certainly cannot expect that the visibihty of the 
quantum effects increases with distance in a naive way. For this reason the resuhs of the 
present paper have no direct relevance to the scenario of modified galactic dynamics in 
EH 120]. 
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A Threshold functions for various cutoff schemes 

In this appendix we collect the relevant properties of the threshold functions $ and $ for 
the cutoff schemes used in the main part of the paper. We employ different cutoff schemes 
in order to gain insight about whether properties of our solutions are independent of the 
scheme used to suppress the IR-modes in the path integral when deriving the beta-functions 
flH.2a| b). This cutoff-scheme dependence is contained in the threshold functions 



mw) =-— / dzz 



w 



(defined for p = 1,2,3,- ■■ and n > 0) through the choice of the dimensionless cutoff 
function R^^\z). Concretely, we employed three different type A cutoffs jTHl E], the 
exponential cutoff, the optimized cutoff [IH] , and the sharp cutoff ^I] . The latter two have 
the advantage that the integrals appearing in ()A.1|) can be evaluated analytically. This 
provides a considerable simplification when solving the RG-equations as for these cases the 
RG-flow of G{k),A{k) is governed by simple first order differential equations. 

A.l The optimized cutoff 

For the optimized cutoff the dimensionless cutoff function R^^\z) is given by 

i?(o)(z)°pt = (1-2)9(1-2). (A.2) 
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Substituting this expression into ()A.1|) and carrying out the integrals we obtain 



These threshold functions are related through the relations 



(A.4) 



A. 2 The exponential cutoff 

The exponential cutoff is given by the following one-parameter family of cutoff functions 

[m Ellin] 

i?(o)(^.,)Exp^ p—- s>0, (A.5) 

exp(s2;) — 1 

where s is a shape parameter. In this case the integrals (jA.lj) cannot be evaluated analyti- 
cally. Therefore we have to rely on a numerical integration when solving the RG equations. 
In order to have a reasonable convergence in the threshold functions we choose the shape 
parameter s = 2. 

A. 3 The sharp cutoff 

The sharp cutoff is defined via the cutoff function 

= = i?0(l-2), (A.6) 

where the limit i? ^ oo is taken after substituting this function into the integrals defining 
the threshold functions. The resulting threshold functions have been determined in 

^\{wf^ = - + w) + <fn, for p = 1 , 

T{n) 

r{n) p-1 {l + w)f ^ ^^7^ 

|.i(u;)'^^=0, forp>l. 
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The <y9„'s are a priori undetermined constants of integration. We choose them as JT] 

ip^ = $;,(0)'^^P(^=i) = n C(n + 1) , (A.8) 



where ({n) denotes the Riemann (^-function. 



B H as a function of the RG scale 

In this appendix we rewrite the Hubble parameter as a function of k, G{k), and A[k). Our 
starting point is eq. ()2.1a|) which (for K = 0) yields 



H^- = ±— [A + 8nGpf\ 



:b.i) 



with the plus (minus) sign corresponding to an expanding (contracting) universe. We then 
substitute eq. ()2.4b|) to obtain H in terms of G, A and their fc-derivatives 



H = ± 



.1/2 



3G' 



- (G"A-GA') 



The k derivatives are conveniently expressed through the modified beta-functions 

~ _l ~ _ 1 

With these new functions the derivatives G', A' are simply 



(B.2) 



(B.3) 



G" = /3G(fc,G,A) , A' = /3A(fc,G,A) 



(B.4) 



Substituting the sharp cutoff ()A.7|) into these equations we find 
kG^ 



where 



Pg 



1 - Bf k^ G 



-51n(l-§)+2C(3)-^r,^ 



= — 

B 

n 



2 " 67r' 



1-Bfk^G' 



51n|l-^)+C(2)-- 



2A 
A:2 



-6 



:b.5) 



:b.6) 
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Eq. ()B.2|) can then be evaluated along any given RG-trajectory k i— > {G{k), A{k)). Thus, 
in terms of ()B.5jl with ()B.6jl . the Hubble parameter for an expanding universe is given by 

1/2 



H{k,G,A) 



A(3G{k,G,A)-G(3^{k,G,A) 
3/3g(A;,G,A) 



:b.7) 
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